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Effective Field Theory for Ultrasoft Momenta in NRQCD and NRQED 
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We propose an effective field theory for heavy quark-antiquark bound states once the soft gluons have been 
integrated out. We also give new results for the matching between QCD (QED) and NRQCD (NRQED). 



1. INTRODUCTION 

Two particle non-relativistic bound states in 
QCD (QED) have, at least, three well separated 
scales: the mass m (hard scale), the typical rel- 
ative momentum \p\ (soft scale) and the typical 
bound state energy E (ultrasoft scale). This al- 
lows to introduce a hierarchy of effective field the- 
ories when sequentially integrating out each of 
these scales. After integrating out the hard scale 
TO QCD (QED) becomes a non-relativistic the- 
ory, the so called NRQCD (NRQED) §. This 
effective theory is local in space and time and it 
is still naturally written in terms of quark fields, 
fntegrating out the soft scale leads to a so far elu- 
sive effective field theory which we call potential 
NRQCD (pNRQCD). We claim that this effec- 
tive theory is local in time but it is non-local in 
space and it is naturally written in terms of wave 
function fields. 

We present some original results for the above 
mentioned effective field theories. In section 2 we 
give the matching coefficients for the four-quark 
operators at one loop for NRQED and also for 
NRQCD in the case of different quark masses. In 
section 3 we put forward our proposal for pN- 
RQCD and give some results for pNRQED. 

2. FROM QCD TO NRQCD 

The matching for NRQCD (and HQET) has 
been known at tree level since long. This can be 



obtained doing the matching at tree level of S- 
matrix elements or just performing a F-W trans- 
formation of the QCD lagrangian. Some results 
at one loop have also been known for HQET. 
Nevertheless, attempts to perform the matching 
beyond tree level using dimensional regularisa- 
tion (DR) in NRQCD have not appeared until 
recently. The problem was that in NRQCD un- 
like in HQET the kinetic term was incorporated 
in the propagator. In DR the high modes are not 
explicitly suppressed by a cut-off /i {p << m) and 
give non- vanishing contributions which break the 
power counting rules. 

Lately several people have addressed this prob- 
lem |2-^] and recently the situation has been clar- 
ified |f|. There, it is claimed that the matching 
should be performed just like in HQET. Let us 
make some comments in favor of this approach. 
The key point is that when doing the matching it 
is not so important to know the power counting 
in the effective theory as to know that the scales 
of the effective theories are much lower than the 
mass. The power counting will tell us the rela- 
tive importance between different operators but 
this would not change the value of the matching 
coefficients. That is, we only need 



TO >> \p\, E, Aqcd 



(1) 



*A.P. acknowledges a fellowship from Generalitat of 
Catalunya. Financial support from CICYT, contract 
AEN95-0590 and financial support from CIRIT, contract 
GRQ93-1047 is also acknowledged. 



whatever the relation between \p\, E and Aqcd 
is. 

In ref. Q DR was used for both ultraviolet 
(UV) and infrared (IR) divergences in the full and 
the effective theory. In fact, it is not so important 
to know the way the UV divergences of the full 
theory are regulated since the comparison is done 



2 



between S-matrix elements which are UV finite. 
Nevertheless, it is essential to regulate in the same 
way the IR divergences of both theories in order 
to cancel. This will happen since both theories 
have by construction the same IR behavior. It 
is also very important, from a practical point of 
view, to regulate the UV divergences of the effec- 
tive theory using DR. In this way, in ref. Q the 
calculation in the effective theory becomes triv- 
ial, since there is no dimensionful parameter in 
the integrand, and the matching coefficients for 
the bilinear terms in fermions are calculated at 
one loop. Nevertheless, the four-quark operators 
were not taken into account. The way to deal 
with these operators is not obvious since we are 
faced with the computation of S-matrix elements 
of four-quarks in QCD and HQET. It is in these 
S-matrix elements, which are never calculated in 
traditional applications of HQET, where the dis- 
tinct IR behavior of two heavy quark bound states 
becomes apparent. Power-like IR divergences ap- 
pear in loops where a quark and an antiquark in 
HQET interact through a potential. We call these 
divergences Coulomb pole. They are naturally 
regulated once the kinetic term is introduced. 

This IR behavior should appear in both the full 
and the effective theory. However, it is important 
to bear in mind that the matching coefficients are 
independent of the IR behavior. Therefore, we do 
not need to regulate the IR singularity with the 
introduction of the kinetic term and hence we can 
take advantage of a more convenient regularisa- 
tion. For instance, a regularisation such that we 
could avoid the effort of computing this pole in 
both theories (which can be very painful). 

The procedure consist of computing the ma- 
trix elements on-shell and at threshold {\p\ = 0). 
In this way there is no scale in the effective the- 
ory, and hence the diagrams in the effective the- 
ory are just trivially zero. Therefore, only the 
diagrams in QCD need to be computed in this 
peculiar kinematical situation. This precise kine- 
matical situation produces IR divergences which 
get canceled by those of the effective theory. The 
computation is done in the MS scheme. 

We stress that the Coulomb pole does not ap- 
pear at all doing the matching in this way. Let 
us explain what happens. In order to define some 



integrals we have to move to dimensions high 
enough for the IR Coulomb singularity to be reg- 
ulated. When coming back to four dimensions 
we can trace back the IR Coulomb singularity as 
poles in dimensions different from four. The point 
is that we have not introduce the relative momen- 
tum and hence DR has no way to reproduce the 
Coulomb pole or any non-local behavior in the 
relative momentum. This fact has already been 
observed in ref. where it was noticed that in 
HQET with a quark and antiquark moving ex- 
actly at the same velocity no imaginary anoma- 
lous dimensions occur using DR. 

The important thing doing the matching is to 
take into account all the non-analytical behavior 
which cannot be obtained in the effective theory. 
We are taking into account all the non-analytical 
behavior due to the masses. The remaining non- 
analytical behavior is encoded in the effective the- 
ory. 

Let us now give the results for the four-quark 
effective lagrangian (non-equal mass case) 
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The QED coefficients are easily obtained from 
these results. We just have to omit dvs and dw 
and replace ^"^^2^ by 1. 

In the equal mass case annihilation processes 
are allowed and they should be taken into ac- 
count. For QED, joining all the contributions we 
get 
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These results are compatible with those found by 
Labelle et al. except for a finite piece in (|^. 

A more detailed explanation of the procedure 
and the full results at one loop for the equal and 
non-equal mass case in QCD will be given in ref. 
§. 

3. FROM NRQCD TO pNRQCD 

In the last section we have integrated out the 
hard gluons. Here, wc integrate out soft gluons. 



with energies of the order of the relative momen- 
tum. 

Two point functions are insensitive to the rel- 
ative momentum and hence the bilinear terms in 
fermions in the NRQCD lagrangian and in the 
pNRQCD at quark level lagrangian will read ex- 
actly the same. However one has to keep in 
mind that in the latter only gluons with ultra- 
soft momenta are kept. On the contrary, four 
point functions do know about relative momen- 
tum and generate non-trivial terms in the pN- 
RQCD lagrangian. These terms are nothing but 
the potential piece of the new lagrangian. Due to 
the massless nature of the gluons the coefficients 
are going to be non-local in the relative space co- 
ordinate although local in time. The important 
point to be realized is that the appearance of a 
potential can be understood as the effect of inte- 
grating out soft gluons. Hence the potential can 
be calculated by matching NRQCD to pNRQCD. 

We have carried out the matching to a given 
order in l/?7i and a using HQET propagators 
and the Coulomb gauge. This produces a very 
strong simplification since the kinetic term can be 
treated perturbatively when computing the po- 
tential. Now, it is very easy to know how far we 
must go in the computation of the potential if we 
want to compute, for instance, the energy up to 
order moi^ . In this case, we must compute the 
matching up to order (^)'*a'', with s, r such that 
s + r < n — 1. The lowest order just gives the 
standard Schrodinger equation. 

In addition we use DR for both UV and IR di- 
vergences, and hence any loop in pNRQCD gives 
zero because there is no scale. The point is that 
on the one hand there is no way of reproducing 
the Coulomb pole since the computation is done 
with HQET propagators in both theories. On the 
other hand any loop with ultrasoft gluons in pN- 
RQCD is zero since these are only sensible to the 
incoming energy and total momentum of the S- 
matrix element, which we set to zero (we remark 
that on-shell quarks in HQET have zero energy). 
Therefore, the calculation from NRQCD leads di- 
rectly to the potential. Further terms in the pN- 
RQCD lagrangian are obtained when matching 
four point functions with an arbitrary number of 
ultrasoft gluon legs. 
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In order to make explicit the distinction be- 
tween soft and uftrasoft gluons it is most con- 
venient to project NRQCD or pNRQCD to the 
two particle sector and promote the wave func- 
tion tlj{xi,X2,t) to a field {ijj{xi,X2,t) is a 3 x 3 
matrix in color space and a 2 x 2 matrix in spin 
space). Then the relative coordinate x = Xi — X2, 
whose typical size is the inverse of the soft scale, 
is explicit and can be consider a small scale com- 
pared with the typical wavelength of the ultrasoft 
gluon which is also of the order of the inverse of 
the energy E. Gluon fields appear in this for- 
malism at the points xi and X2- Ultrasoft glu- 
ons are those which B^{xi, t) can be expanded 
about the center of mass coordinate in a power 
series of derivatives and relative coordinates, the 
so-called multipole expansion (refs. |^,^ also deal 
with ultrasoft momenta through the multipole ex- 
pansion). Then the most natural way of writ- 
ing the pNRQCD lagrangian is as a functional of 
tjj{x, X, t) and B^{X , t), which is local in Vj, Vjj 
but non-local in x. For definiteness we will give 
LpNRQCD at the leading order in the potential. 



In this way (H) reads 



^pNRQCD = 



(fx(fXdttryTp\xi,X2,t) (9) 
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\Xl-X2\ \ 

Recall that the gluon fields in the covariant 
derivatives are ultrasoft and hence they must 
be multipole expanded which spoils the explicit 
gauge invariance. However, this can be recov- 
ered by writing 'ijj{x,X,t) in terms of the scalar 
S{x,X,t) and octet 0(x,X,t) wave function 
fields, whose local gauge transformations depend 
on the center of mass coordinate only. Namely, 

iP{xi,X2,t) = P[e^^2 ]S[x,X,t) 
+P p Q ^'^--j 0(f , X, t)P[e' ^'''] 
ip{xi,X2,t) g{xi,t)ip{xi,X2,t)g^'^{x2,t) 
S{x,X,t) S{x,X,t) 
0{x, X, t) ^ g(X, t)0{x, X, t)g-\X, t) 
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This lagrangian suffices to obtain the leading 
non-perturbative contributions to the two heavy 
quark bound states when E >> Aqcd @| (see 
also @). 

If we leave aside non-perturbative effects (~ 
Aqcd), each term in the lagrangian above has 
a well defined size, unlike in the NRQCD la- 
grangian. Relative coordinates x and its asso- 
ciated momentum must be counted as soft 
scales {{x)~^ ~ ^ ma). Gluon fields B^{X,t) 
and derivatives with respect to the center of 
mass coordinate Vjj must be counted as ultra- 
soft scales mc?). Then if we wish to cal- 
culate a given observable to a given order in a 
we know immediately which terms are to be kept 
in the pNRQCD lagrangian. As an example, let 
us present LpjsiRQED from which one can obtain 
next to next to leading corrections to the energy 
(~ TOa^). 

CpNRQED= j (f'xdPXdtS\x,X,t) 
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The coefficients q are given in [Q and dss, dsv in 
(0) and (|). c, = l+0{a) and = 0{a)+0{a^) 
are obtained from tlie one loop matching between 
QED and NRQED. The last term corresponds 
to the ultrasoft photons which contribute at this 
order. The potential terms are obtained upon 
matching NRQED to pNRQED up to one loop. 
6V encodes the computation of the potential at 
one loop, ft reads 

SVix) = (s^^^x)lni^^ + 7?-reg^ 
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with unknown C. The first In v (which was pro- 
duced by an UV divergence in the potential) gets 
canceled by the In v in dss ■ The second In i/ (with 
IR origin) cancels with a piece of the UV diver- 
gent contribution coming from the ultrasoft pho- 
tons. The remaining contribution from the ul- 
trasoft photons cancels the contribution coming 
from C]j . The net result is the total energy to be 
scale independent as it must be. 

With ( pT| ) the binding energy at order ma^ for 
arbitrary n, I states can be calculated. We find 
agreement with for the hyperfine splittings. 
The ma^ In a correction has also been calculated 
finding agreement with known results pl] |. We 
have also used these techniques with DR to re- 
produce the Lamb shift in the simpler case of a 
hydrogen-like atom [T^ . 

We remark that in this section we are assuming 
\p\ » E which was not needed in the previous 
section. As far as we are making the matching 
to some order in a (we are computing the poten- 
tial perturbatively) we also assume |p| >> ^qcd- 
Notice that the relative size between E and Aqcd 
is left arbitrary. We can further distinguish be- 
tween two situations: (i) E >> Aqcd and (ii) 
Aqcd E. For the situation (i), as we men- 
tion before, we can calculate perturbatively from 



pNRQCD and parametrise the non-perturbative 
contributions by means of local condensates. For 
the situation (ii) the calculations in pNRQCD 
cannot be carried out perturbatively anymore. 

If Aqcd > IpI the matching to pNRQCD can- 
not be carried out perturbatively. Even in this sit- 
uation NRQCD is extremely useful to parametrise 
non-perturbative contributions in many processes 

§• 
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